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We theoretically investigate the Carlson-Goldman (CG) mode in two-dimensional clean d-wave 
superconductors using the effective "phase only" action formalism. In conventional s-wave super- 
conductors, it is known that the CG mode is observed as a peak in the structure factor of the 
pair susceptibility 5(f2, K) only just below the transition temperature Tc and only in dirty sys- 
tems. On the other hand, our analytical results support the statement by Y. Ohashi and S. Takada, 
Phys. Rev. B 62, 5971 (2000) that in d-wave superconductors the CG mode can exist in clean sys- 
tems down to the much lower temperatures, T « O.lTc. We also consider the manifestations of the 
CG mode in the density-density and current-current correlators and discuss the gauge independence 
of the obtained results. 

PACS numbers: 74.40.-(-k, 74.72.-h, ll.lO.Wx 



I. INTRODUCTION 



More than 25 years ago, an unusual propagating sound-like {fl/K ~ vqq) Carlson-Goldman (CG) mode in charged 
superconducting systems with the velocity vcg = 10'^ — lO^m/c was discovered (see also [Q). It was a widespread 
opinion before the CG mode discovery that since in the charged systems the sound-like Bogolyubov- Anderson (BA) 
^ mode associated with neutral superconductors is converted to the plasma mode due to the Anderson-Higgs 
mechanism, there is no sound-like phase mode in charged systems. 

A magnificent effort (see, for example, the reviews ^, 0, ^) was made to understand the mechanism, responsible 
for the appearance of the CG mode and its relation to other phenomena of non-equilibrium superconductivity. While 
the majority theories of the CG mode [|| ^,0, ^) were essentially based on the kinetic equations which are usually 
derived using the quasi-classical Green's functions, the paper by Kulik, Entin-Wohlman, Orbach ^ used a more 
conventional approach based on the Matsubara Green's functions without kinetic equations. In the subsequent papers 
of Takada with coauthors (see ||l^, |ll| and references therein) the approach of ||] was further developed and very 
recently applied for the case of d-wave superconductivity . The collective oscillations in d-wave superconductors 
were also studied using the kinetic equations for Green's functions (see also |Q). 

Since the discovery of high-Tc compounds, the dj;2-y2 superconductivity has attracted much attention ||l5|| and the 
claim of ||l^ that the CG mode in clean d-wave superconductors may survive in a much wider region of temperatures 
down to 0.2Tc appears to be very different from the established properties of the CG mode in s-wave superconductors, 
so that it is important to check it by an independent calculation. 

On the other hand, the importance of phase fluctuations in high-temperature superconductors (HTSC) (see, e.g. the 
review [p^), stimulated interest in the derivation of the "phase only" effective actions from the microscopic theory. It 
is important to emphasize that although there is no commonly accepted theory of HTSC, it seems reasonable that one 
can use a simple BCS-like approach to describe the properties of HTSC below the critical temperature, Tc even though 
such an approach fails above Tc- Relying on this argument the "phase only" actions for d-wave superconductors were 
recently derived in |l8[ However, the only phase excitations which are described by these actions are the BA 
mode in the neutral superconductor p8[ | and the plasma mode ||l^, which appears when the Coulomb interaction 
is taken into account. This corresponds to the standard paradigm which does not yield the existence of the CG mode. 
Thus the purpose of the present paper is to investigate which ingredient is missing in the treatments of |ll7| , p^ , , 
so that the CG mode does not appear in these approaches and to establish a link between the results of |^ and the 
"phase only" action formalism. This missing link is established here and the CG mode is obtained within the effective 
action formalism. Our main results can be summarized as follows. 

1. We extend the "phase only" effective action formalism for charged systems to incorporate the density- current 
coupling which was so far considered only using other methods [|, [l^, y, |l| (sec also [|o[ 1^, |2|, where the 
effect of this coupling appears to be important for the description of dirty superconductors). In neutral systems the 
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density-current correlator does contribute in so called Landau terms of the effective action [ [18| , p4[ , so that the correct 
expression for these terms can only be obtained when this correlator is taken into account. 

2. We show that when the density-current coupling is included it becomes possible to obtain the CG mode using 
the "phase only" action. In particular, we derive an analytical expression for the pair susceptibility structure factor 
and solve numerically the equation for the CG mode velocity. 

3. We show the gauge independent character of the equation for the collective phase excitations one of the solutions 
of which is the CG mode. Establishing a link between the pair susceptibility and this gauge independent equation for 
the phase excitations we argue that the peaks in the structure factor associated with the CG mode are independent 
of the choice of the gauge. The gauge independence of the equation for the CG mode used in the previous papers 
pO| , |ll| , is also shown applying the identity derived recently in [ p2| . 

4. We consider possible manifestations of the CG mode in the gauge independent density-density and current- 
current response functions. 

5. We derive analytical expressions for the density-density, current-current and density-current polarization func- 
tions for 2D clean d-wave superconductor at T ^ Tc. 

The paper is organized as follows: In Sec. II we introduce our model. In Sec. Ill, we describe all details about the 
formalism used in the paper, necessary for further understanding. The various forms of the effective actions (for the 
phase and electric potential, "phase only" and "electric potential only") are expressed in terms of density-density, 
current-current and density-current polarization functions in Sec. IV. The general expressions for these polarization 
functions are given in Appendix their derivation for d-wave superconductors is considered in Appendix ^ and the 
nodal approximation employed for this derivation is briefl y d iscussed in Appendix In Sec. |l^ we also discuss in 
detail the difference between the present and other ^ |l^, ]l^, |2[ |2l|, |2[ approaches. In Sec. we briefly recall 
the properties of the phase excitations in the absence of the density-current coupling and stress some points like 
gauge independence of the equation for the collective phase excitations and the properties of the gauge independent 
density-density and current-current correlators which are pa rticul arly useful for better understanding of our main 
results which are presented in Sec. VI. In particular, in Sec. VIA we derive the equation for the CG mode, give its 
physical interpretation and discuss the gauge independence of th e present and previous approaches. In Sec. VI B we 
present the results for the velocity of the CG mode and Sec. VI C is devot ed t o the structure factor (the calculational 
detail for these sections are given in Appendix ^). We conclude in Sec. VII with a discussion and summary of our 
results. 



II. MODEL HAMILTONIAN 



Let us consider the following action (in our notations the functional integral is expressed via e'^) 



S = - I dT V / d\iji{T,v){dr - ieAo{T,r))i;,{T,r) + H{t) 
Jo I „ J 

where the Hamiltonian H(t) is 

^(^) = E / rf''^V'i(T,r)[£(-zV-^A(T,r))-^]V'.(T,r) 

+ \j d\, j d2r2(^^t(^^ri)V.(r,ri)-n)14(ri-r2)(X]Vl,(T,r2)iA..(T,r2)-n) 



(1) 



(2) 



Here ipa-iT, r) is a fermion field with the spin a =1 , |, ct = —a, t is the imaginary time and V{vi] V2) is an attractive 
short-range potential, 14 (ri — ^2) is the long range Coulomb interaction, n is the neutralizing background charge 
density. Throughout the paper we call the superconducting system neutral if the last term of Eq. ^ is omitted and 
charged if this term is taken into account. Even in the latter case the whole superconductor remains, of course, neutral 
due to the neutralizing ionic background. 

We assume that the momentum representation of V{yi\Y2) contains attraction only in the d-wave channel (see 
the discussion in [0). The Fourier transform of the Coulomb interaction 14 (q) depends on the detail of the model. 
It can, for example, be taken T4(q) — Aire^ /q^ in 3D, T4;(q) = 27re^/g in 2D or more complicated expression (see 
e.g. |l^) if the layered structure of HTSC is taken into account. However, the detailed expression is not crucial 
for the CG mode, because the mode appears when the Coulomb interaction is screened out by the quasiparticles. 
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(The form of the expressfon would be, of course, essential for the analysis of the plasma mode |]T^, ^.) The form 
of dispersion law, e(k), is also not essential because the final results for the rf-wave case will be formulated in terms 
of the non-interacting Fermi velocity -vp = 9£(k)/9k|k=kF ^nd the gap velocity va = 3A(k)/9k|k=kF: where A(k) 
is the momentum dependent superconducting gap. We will also use the parameter ao = vf/va which is called the 
anisotropy of the Dirac cone. Throughout the paper H — ks ~ 1 units are chosen. An external electromagnetic 
field A — (Aq, A) was introduced in the action Eq. (|^) to calculate various correlation functions using the functional 
derivatives with respect to this external source field. 

III. DERIVATION OF THE EFFECTIVE ACTION AND THE STRUCTURE FACTOR 

The derivation of the effective "phase only" action for neutral (see e.g. |l^, and charged jl^, |l^, ^ s- 
and d-wave superconducting systems is widely discussed in the literature, so we briefly recap the main steps, including 



the functional integral representation for the structure factor, and making in Sec. IV a point on the appearance of 
the term which couples density and current. 

The first step of the derivation is to use the appropriate Hubbard- Stratonovich transformations to decouple four- 
fermion interaction terms in the attractive and repulsive channels. The attractive part of the interaction was recently 
considered in detail in Sec. H of using the bilocal Hubbard-Stratonovich fields $(t, ri;r2) and <I>''^(t, ri; (see 
p7[ for a review), so we show explicitly the corresponding transformation only for the Coulomb interaction: 



exp 



2 „ „ 



I [ dr [ dVi / dV2(V'i(r,ri)V'a(T,ri) -7i)K(ri,r2)(V';t'(^'i"2)?A<,'(^,r2) -n) 

^eif{T, ri)V;"^(ri,r2)e(y9(T, r2) - ie(^(T, ri)(V'i(T, ri)'(/'a(T, ri) - n)S{ri - r2) 

where the Hubbard-Stratonovich field has the meaning of the electric potential. Thus, the partition function is 
Z=J V^''V^V<^''V<i>Vtpexp[S{^\^,^'',^,ip,A)], 



(3) 
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S{^1^,<P\^,^,A)^ dr J dVi J dV2|-:^^^^^|$(r,ri;r2)|2-^e^(T,ri)y-i(ri,r2)e^(T,r2) 

+ *t(r,ri)[-9^ + icTsAo + ieraipir, ri) - Ta^i-iV - r3-A)]*(r, r2)(5(ri - ra) 

c 

-ie(p{T,ri)nS{ri ~ ra) + $^(7-, ri;r2)*^(r, ri)T_*(r, v-^) + *^(r, ri)r+*(T, r2)$(r, ri;r2)} 

where ^ and are the Nambu spinors, $(— irsV) = e(— iraV) — /i and ra, t± — (t\ ± iT2)/2 are Pauli matrices. 

To consider the Hubbard-Stratonovich field $ it is convenient to use the relative r = ri — r2 and center of mass 
coordinates R = (ri -I- r2)/2, so that $(r, ri, Y2) = $(t, R, r). Now, using the functional integral representation, the 
imaginary time pair susceptibility is defined as 

7;(r,Ri -R2) = - / I?*tp^2?$tx»$X)^$(r,Ri,0)$t(o,R2,0)exp[S'(*^*,$t^$,(^)]. (5) 

Z[A — \)) J 

Since the distance |Ri — R2I is expected to be larger than the internal Cooper pair scale, it is possible to put r = 
in Eq. (H). The structure factor 5(57, K) which used to present the experimental data ||, |j is related to the real 
frequency pair susceptibility ry(i51„ — s- 57 iO, K) by 

S{n, K) = -2[1 - cxp{-/3n)]-^lmr]{n + iO, K) w -—lmT]{n + iO, K) . (6) 

The definition of the pair susceptibility Eq. (||) is apparently gauge dependent, since the auxiliary Hubbard- 
Stratonovich field is gauge dependent. Nevertheless, as we discuss later the poles of r]{n, K) are gauge independent 
and this justifies the use of Eqs. (I) and (|) to extract the observable values. 

The simplest way to study the low energy phase dynamics is to employ the canonical gauge transformation 

/ ie(x)/2 Q \ 

*(^)- f „ g-,.(.)/2)T(x), x^(r,r) (7) 
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separating the phase of the ordering field 

$(t, R, r) « A(t, R, r) exp[ie{T, R)] . 
Then after tlic integration over the Fermi-fields the partition function becomes 



Z = y AVAV9Vipcxp[-f3n{A,de,ip, A}] 



(8) 



(9) 



where the effective potential 

pn{A,de,^,A} = 



A2(T,R,r) 1 



^(ri-ra) 2 



+ -e(^(T,ri)T/ ^(ri,r2)e(/7(T,r2) + ieip{T,ri)nS{ri - r2) 



with the inverse Green's function 



TrLnG- 



.1 (10) 



(11) 



(Ti,ri|S(a0,^,A)|T2,r2) 



6 

T3 ( ieip{Ti, ri) - ieAo{Ti,ri) ) +T30i(ri) +/02(ri 



'^(ti - T2)(5(ri - r2) . 



(12) 



For T <C Tc it is reasonable to neglect the amplitude fluctuations and assume that the amplitude of the order 
parameter A(t, R, r) does not depend on R. Then the frequency-momentum representation of G in Eq. ( pl| ) is the 
usual Nambu-Gor'kov Green's function 



t/(iw„,k) 



(13) 



where, because d-wave pairing is considered A(k) = A/2(cosfc2:a — coskya) (a is the lattice constant) and a;„ — 
7r(2n -I- 1)T is fermionic (odd) Matsubara frequency. Since in what follows only the low temperatures, T ^ A(r) are 
considered, we can replace the temperature dependent amplitude A(r) by its zero temperature value, Aq = A(T = 0). 
Thus all linear low temperature dependences of the polarization functions considered below are due to the nodes of 
A(k), but not to the temperature dependence of A(r) itself. 

The precise form of the operators Oi and O2 in Eq. ( [1^ ) which depends on the particular form of the tight-binding 
spectrum £(k) is given in jl^ (see also |19 for the formulation of the general rules for representation of S). It is 
essential, however, that the coordinate representation of S does not depend on the phase 9 itself and contains only 
its derivatives. Thus the coordinate representation of ilkin is also expressed via the derivatives of 9. This property 
is particularly convenient for studying 2D models when a constant space independent phase is prohibited by the 
Coleman-Mermin-Wagner-Hohenberg (CMWH) theorem. 

Since we are interested only in the phase dynamics in the presence of Coulomb interaction, in what follows we 
consider only the phase 9 and the electric potential ip dependent parts of the thermodynamical potential Eq. (ITT 
This part of f2 which we denote as Qkin can be present as a series 



nkin{d9,ip} = dr 
Jo 



d^n / 



-etp{T, ri)V^ ^(ri,r2)ev3(r, r2) + ieip{T,ri)nS{ri - r2) 



00 

■E 

11=1 



TTr > ■ -(^E)" 



(14) 



This way of deriving the effective action has many advantages. The main among them is that the gauge invariant 
combinations hdT9/2 — eip— bAq and hV9/2 — e/cA. are explicitly present during all stages of the derivation |Q. This 
property is obviously related to the introduction of the phase via the gauge transformation Eq. (^. There is no need 
because of this to keep the external electromagnetic field ( , A) during the intermediate stages of the derivation 
since it can be easily restored following the above mentioned prescription which in the frequency-momentum space 
are 



inn9{K) -> in„9{K) + 2eAa{K); iKa9{K) 



iK^9{K)--A^{K). 
c 



(15) 
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Differentiating with respect to this source field we will derive physical correlation functions (see the discussion in [pT|) 
in what follows. It has to be stressed that the minimal coupling prescription ([l5| ) does not guaranty itself the gauge 
independence of the final result. The gauge independent treatment of the transformations (0) and (||) for d-wave 
pairing is, in particular, one of the complications p9|. 

In the previous studies of the CG mode [lO|, |ll|, 12 1 the phase field was introduced using the expansion of the 
ordering field ^{x) around the equilibrium value A via ^{x) = A + ^i{x) + i^2{x) and associating the fields ^i{x) 
with the amplitude and ^2{x) (or to be more precise $2(a;)/Ao) with the phase fluctuations. Although, as will be 
discussed below, the final result obtained in the both methods agrees, the present method of the investigation of the 
CG mode is more transparent because it explicitly uses the gauge independent combinations of the fields over the 
whole derivation. For example, one can easily recognize that S contains a gauge independent Cooper pair chemical 
potential 6/2 — e^p [|| which in other approaches has to be collected from the different parts of the equations. Other 
advantages of the present approach will be discussed in the subsequent sections where the effective action is presented. 

To consider the phase and charge fluctuations at the Gaussian level it is sufficient to include only the terms with 
n = 1 , 2 in the infinite series in Eq. (H). Finally, we rewrite the pair susceptibility Eq. (||) in the new variables 

?7(r,Ri -K2) = - J AVAV9VifA{T,Ri,O)exp[ie{T,Ri)]A{0,R2,O)exp[-ie{T,R2)]exp[-f3n{A,de,if}] . (16) 

As was mentioned above, we are interested only in the phase fluctuations structure factor neglecting the presence 
of the amplitude fluctuations. This implies that one can use a saddle point approximation for A, so that omitting 
unimportant constant in the expansion of the exponent in Eq. (|l^) one arrives at 

r,(r, Ri _ R2) = J V0Vip0iT, Ri)0{t, R2) exp[-/3^2fe,„{a0, ip}] . (17) 



Writing Eq. (|17|) we also expanded the exponents which were present in Eq. (^^, because there are no free vortices 
in the system for T < Tc and the multivalued character of the phase is irrelevant. This approximated form of the 
pair susceptibility is equivalent to the expressions for the susceptibility used in Refs. j9| |l^, |ll| |l^. Expanding the 
exponents we neglect a widening of the structure factor peaks which is related to the absence of the long-range order 
in 2D (CMWH theorem). It is known, for example, from the analysis of the dynamic structure factor of lattices 
pl| , ^ that (5-function phonon resonance obtained in 3D harmonic crystals in 2D case is converted to the power law 
singularity 

where /(a) is a function of a which goes to zero as a — > 0, so that in this limit the structure factor transforms to 
(5-function. Since for low temperatures a(T) <C 1 we may safely neglect this effect of widening because it does not 
move the position of the peak and we are primarily interested in the temperatures T <C Tc- 

IV. GENERAL FORM OF THE EFFECTIVE ACTION 

In this section we present the effective potential r2fei„{0, 1^} and discuss the term which leads to the appearance 
of the CG mode. We also derive the effective "phase only" and "electric potential only" actions integrating out the 
electric potential ip and the phase 6*, respectively. 

A. The effective action and polarization functions 

Calculating the terms with n = 1, 2 in Eq. ( p^ ) (see e.g. ^) one arrives at 

(3nk^n{e, - g ^ J ['pi-KjAe^V.'^KMK) + (^^}„0(-if ) - 2e^{-K))Il:i3{K){in„e{K) + 2e^{K)) 

+ e{-K)[Kf + ni^{K)]K^Kpe{K) 

+0(-K){-K^)U^oiK){i^n0{K) + 2ev{K)) + {i^l^e{-K) - 2e^{^K)){-K^)U'^^{K)e{K)] , 

(19) 



n— — 00 
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where we introduced short-hand notations K = (iri„,K) with K being 2D vector (summation over dummy indices 
a, /3 = 1,2 is imphed). In Eq. (^9|) the current-current polarization function, IIoo is 



oo „ 



1— — 00 ' 



(2^ 



7roo(jf^„, K; iw;, k)u_Fa(k)wF/3(k) 



with the Fermi velocity ^^^^(k) = dS^i\<i) / dk^', the density-density polarization function, 1133 is 



°° [ (fk 

n33(if]„,K) = T ^ / — — 2 7r33(if2„,K;iw/,k 

— 00 



^03 

oc 



n^3(zfi„,K) = T ^ y ^7ro3(zr!n,K;ztu;,k>^^„(k) 



and the density-current polarization function, Ilgg is 

TTy in Eqs. (|o|) - (^) is given by 

TTy (iSl„, K; icj;, k) ee tr[t/(iw/ -|- k + K/2)TjfJ(zcj;, k - K/2)rj] , 

and Aq^ in Eq. (|l9|) is the first order contribution in the superfluid stiffness: 

d^k 



(2n) 



i(k)m;i(k). 



-\k) = d^a^)/dk^dk^ 



with 



n(k).l-mtanh.^('^) 



E{k) 



2T 



Eik) = v/C'(k)+A2(k) 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



Writing Eq. ( [l9| ) we omitted the linear time derivative term (see e.g. |T^, |T^) which is irrelevant for the present 
analysis. 

The general expressions for the polarizations (^ ) - (|2^) are given in Appendix |^ (see also (l8|) and calculated 
analytically for 2D clean d-wave superconductor in Appendix |^ (a brief discussion of the nodal approximation used 
to calculate these polarizations and the transformation to the global coordinate system are given in Appendices ^ 
and|^, respectively). 

As an example we show in Fig. |l]the real part of the density-density polarization functio n, 11 33 (f^, K) (this function 
is in fact just the Lindha rd's function for the superconducting state) which is given by Eq. (C3) (and its dimensionless 
form by Eqs. (C4) and (C5)) as a function of fl/vpK for the different directions of K. The angle </> determines the 
direction of K with respect to the (lO)-direction in a such way that — tt/A corresponds to the nodal direction (see 
also Eq. (D5) and the explanation in Appendix Comparing this figure with Fig. 7 from ||l^ (our definition of 
is equivalent to the definition of the angle 6>q used in ||l^), which was obtained by numerical integration of Eq. (|Al| ), 
one can see that our analytical expression ( p3[ ) gives essentially the same result. In particular. Fig. |^ shows that 1133 
has a peak at U/vpK — cos{(f) — 7r/4). Furthermore, as shown, for instance, for (j) — 7r/16 there is a lower peak in 
Ren33(fi, K) at n/vpK = cos{4> — 7r/4 -|- 7r/2). Note that the cases = and (j) = 7r/4 are "degenerate" because for 
(j) = the lower peak coincides with the upper one and for = 7r/4 the lower peak is at Vl/vpK = (see also Fig. ||, 
where the case 4> = 0.2, 0.237r is shown). 

In the origin of these peaks was related to the gap nodes which can be regarded as two "one-dimensional normal 
state electronic bands" toward 4> ~ ±7r/4. These "normal bands" are able to screen out the Coulomb interaction in 
certain regions of the Fermi surface even for T <ti Tc and t his s creening along with the presence of nQ3 will make 



the appearance of the CG mode possible. Substituting Eq. (D5) in the analytical expression ( |C3|) it is indeed easy 
to see that these peaks are due to the square root singularities of 1133. These are the same singularities which are 
present in 2D normal state Lindhard's function due to the lowered dimensionality of the momentum integration 
but since d-wave superconducting state is considered the position of these singularities does depend on the direction 
of K with respect to the Fermi surface. Finally we note that these singularities in Ren33 at 17 = vpK cos((/i— 7r/4) + 
is accompanied by the singularity in Imn33 at = vpK cos{(j} — 7r/4) — which was considered in p8[ . 
The effective potential Eq. ( [T9| ) becomes more tractable in the matrix form 



T 



E 



dK 

(2^ 



[e{-K) eip{-K)]M-'^ 



e{K) ■ 

eip{K) 



(26) 
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where 



-2ir!„n33(if) + 2x„nfo(if) Ai-u^^iK) + v-\k)) 

with the bare (unrenormahzed by the phase fluctuations) superfiuid stiffness A"'^ = A^'' + n^l^{K). 



(27) 



B. Comparison with other approaches 



Let us compare our effective action Eqs. (|2^), ( p7| ) with the action obtained in and see the differences between 
the present and previous |ll|, |l^ approaches. As one can notice, the only difference between Eqs. (^), (|2^) and 
Eqs. (25), (26) in is due to the presence of the density-current polarization function, It is a general belief 

that this term has to be zero due to the "symmetry arguments" ]25[| . However, as shown in (see also [|l8|, where 
the d-wave case is considered) this correlator has to be taken into account to obtain the correct expressions for the 
Landau terms of the effective action. This is the term which nontrivially couples phase and density fluctuations and 
makes the CG mode possible in the present approach. 

From this point of view, the role of no3 is the same as the role of the phase-charge coupling 



n23(ia.,K)=r J2 J ■^7r23{tnn,K;iLOi,k)j{k), 



7(k) 



2A(k) 



(28) 



with 7r23 given by Eq. 
1123 = -1132, while 1X03 



(23) in the approach of Kulik et al. |^ and Takada with coauthors ||T^, |Tl|, |T^. Note that 



TTQ 



It is interesting that techniques essentially similar with |9[ |T^, |T^, have been used in |Q and ^ to 
consider suppression of the critical temperature in disordered superconductors. In |2j, |2^ both amplitude and phase 
fluctuations were taken into account and to consider the influence of nonmagnetic impurities the electronic Green's 
functions had 4x4 matrix structure. The main difference between 12, ^ nevertheless, remained 

the same, viz. the order parameter phase was expressed via the operator O2 — ^'^2^', as summarized in Table II of 
pl|. T hus T2 also enters the phase-density correlator, 1123 in the notations of [ij] (or II^p in the notations 

offii,!!). 

In our opinion, the physical meaning of 1123 is, however, more obscure than that of Ilpg. Indeed, since 1123 is 
expressed via the Pauli matrix T2, so that it seems like the phase itself is a dynamical variable on its own, while 
physically meaningful are only the space and time derivatives of the phase. 

These derivatives can only enter into the formal expressions as a current via / and as a density via T3 matrices, 
respectively. This property is obviously present in the definition Eq. ( p^ ) of nQ3 which thus has the more clear meaning 
of a density-current polarization function. 

Another important difference between the present and previous ||, |l^, |ll|, |l2[ ^ |2^, ^ approaches is that the 
present derivation does not need an explicit use of the gap equation for A. For example, in |S| the charge conservation 
follows from the explicit use of the gap equation, while in the present approach as we will discuss later, the charge 
conservation is already built in the phase dynamics itself. 

The fact that the present approach does not rely on the particular form of the gap equation is more convenient 
for modehng HTSC, where the gap does not close at the critical temperature Tc, so that the equation A(T°) — 
gives only the mean-field transition temperature, T^. Thus another definition of the true critical temperature, Tc is 
necessary. As recently discussed in |Q (see also ) , i t is r easonable for HTSC to estimate Tc as the temperature of 
the Berezinskii-Kosterlitz-Thouless transition (see Eq. ( |C1C ) in Appendix |c|) . 

Finally, it is worth to mention here the recent papers |35|, where the approach very similar to that of the present 
paper was employed to study the CG mode in the model of color superconducting quark matter. One of the advantages 
of [^5| is that it treats the problem in an explicitly gauge invariant way, while here, the Coulomb gauge is already 
imposed in writing the Hamiltonian Eq. (^ and then, when necessary, the gauge independence of the results obtained 
is discussed in a more intuitive way. In general, however, to prove the gauge independence of the results an arbitrary 
gauge has to be considered, to show that the physical observables do not depend on the gauge fixing parameter as 
done in (see also |Q , where the gauge invariance of the physical quantities calculated using the T = effective 
potential is discussed). Following this route one can obtain the "relativistic" (see the second paper in for the 
details of the proof) generalization of Eq. ( |l9| ) which contains a gauge fixing parameter A 



pnk^n{o,A^}=Tj2 



1 ( 1 



8 TT 



(29) 



+ {-iK^,e{-K) - 2eA^{-K))n^,{K){tK,0{K) ~ 2eA,{K))} 
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with = (17„, K), = K^K^ = nl + K2, /X, = 0, 1, 2 and Aq{K) = -ip{K) and c = 1. Note that in Eq. (Eg) we 



have the whole electromagnetic potential instead of the Coulomb component present in Eq. (19). The polarization 
tensor Ilfj_^{K) is obviously related to the polarizations used in Eq. (|T^). The question of gauge independence (or 
dependence) can be addressed considering how the calculated quantities depend on A. 

C. Effective actions for the phase and electric potential 

Integrating out ip and 9 from Eq. ( |26| ) one can obtain, respectively 

/3f^fe^„{n-| E j ^0{-K)Mg\in^,Vi)e{K), M,^ ^ M^,^ - M22M^^HI^,^ (30) 

n— — oo 

and 

n— — oo 

It is evident that = e^M'^MuMg'^ = M^^ deiM^^ , so that if M~^,M~^ ^ the equations Mg^ = = 

detAl"^ = are equivalent. 

Using Eqs. ( p7| ) and (|3^) it is straightforward to obtain that the structure factor Eq. (^ is given by 

5(ri,K) = ^^ImAle(17 + zO,K). (32) 

V. PHASE EXCITATIONS IN THE CHARGED SYSTEM IN THE ABSENCE OF THE 

DENSITY-CURRENT COUPLING 

A. Equation for the collective phase excitations 

Let us assume for a moment that there is no density-current coupling (Ilgg = Ilgp = 0) and discuss briefly the 
collective excitations which follow from Eqs. (^^, (|3^) and (|3l|). As mentioned above the matrix M'^^ reduces in 



this case to the known expression |17 . Therefore, it is not surprising that the "phase only" action Eq. ( pO| ) takes the 
form 

M,\in^,K.)^~nl%:, + K'^PK^Kp, ^ —^^^^1-^^ , (33) 

which coincides with the corresponding expression in (see also |2^). The dispersion law of the collective phase 
modes is defined by the equation 

>l^i(ia» ^ + iO,K) = 0. (34) 
This equation can also be regarded as a direct consequence of the charge conservation 

dp{t,v) 



dt 

where the current and charge density are defined via 



V-j(i,r)=0, (35) 



' = ' SA„ ' ^^^^ 

where the electromagnetic field A in flkin was restored using the rule Eq. (|l^). Evaluating Eq. ( ^ ) one arrives at 
Eq. (Bsh with Ilpg — 0, so that Eq. ( ^^ is indeed recovered. There is, in fact, no surprise in the link between Eqs. ( |35| ) 
and ^J) which is just the consequence of the way how we introduced the phase 9 in Eq. (^. 

As discussed in |2^ the only solution of this equation for K ^ is the plasma mode. Using as an example the 3D 
form of the Coulomb potential and assuming that A"^ = A6ap which is valid for isotropic system with niafi — mSafi, 
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one obtains from this equation that the plasma frequency, LOp = v47rAf? for the hmit K ^ 0. The expression for Up 
can be reduced to the standard 47me^/m if one uses the superfluid stiffness K — n/m obtained for the continuum 
model with s-wave pairing. 

It is clear that for the plasma mode Q.p{K.)/K cx3 as K — > 0. This property remains valid even if 2D Coulomb 
potential is used and it makes plasmons different from any sound mode with ^/ K const as K ^ 0. For example, if 
a neutral superconductor were considered the polarization Has in Eq. ( |33| ) would be replaced by Has and the solution 
of Eq. (Is!) for K ^ is the sound- like BA mode and its velocity is given by Eq. (C9). 



B. Phase excitations via electric potential propagator and gauge independent density-density and 

current-current correlators 

It is instructive also to look at the form ^A^p for the electric potential ip which is 

>.,-(.»,.,K,.v( ,^;;;^^^^ .K-(K)), ,37, 

Considering the same example of isotropic system with 3D Coulomb potential Eq. ( p7| ) in the limit K ^ and 
U/K ^ oo can be reduced to the known expression (see e.g. eF 




A^-i(za.,K) = — 1 + . (38) 



It is obvious that the discussed above plasma mode can be also seen in A^;^. Indeed after the analytical continuation 
iU,n ^ f2 + iO, Alip(r2,K) acquires a pole at = uop. 

It is also useful to evaluate gauge independent density-density and current-current correlators which are defined as 

K)- ^'T'-^i^} ^^^m K)- ^'T'-^i^} (39) 

where 

Z{A}^ J VeVipeM-f3^krn{e.ip,A}] (40) 



with the external field (Aq, A) restored using the rule Eq. (15). Then we arrive at the standard expressions |l 



and 

Again assuming that A"'^ = ASa/s one can reduce Eq. ( ^ ) to the known form of the current-current correlator 

-nlu,,s^0 + AK'(s^,-^) 

x"^{inn,K) = e^A ^ ^. (43) 

The difference between gauge independent Eq. (^l]) and Eq. (^) and gauge dependent ftaa, A"'^ in Eq. ( |33| ) was 



recently discussed in MJj. The link considered above between charge conservation Eq. (35) and Eq. ( p4[) shows that 
the solutions of Eq. (Q) are gauge independent, even though Ilaa, A"'' are gauge dependent. At the more formal 
level it can be argued that Eqs. (|4l| ) and ( ^2|) are gauge independent because even if one derives them starting from 
Eq. ( p9| ) the gauge fixing parameter A does not enter the final result and one obtains the same expressions. 

The putting this in another way, one can say that the position of zeros is gauge independent, because 



these zeros coincide with the poles of the gauge independent x(^jK). As we will see in Sec. VI Dj , density-current 



coupling modifies both Eqs. ( |4l| ) and ([42]), nevertheless the general argument about zeros of Mg (^^,K) (poles of 
A^e(^^,K)) remains valid. 
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VI. PHASE EXCITATIONS IN THE PRESENCE OF THE DENSITY-CURRENT COUPLING 



Here we generalize all expressions from the previous section for nonzero 1130 • In particular, Eq. ( |33| ) becomes 



Ti^,{KWUK) (44) 

i-n33(W(K)' +-n33(if) + K-^(K)- 



The dispersion law for all collective phase excitations is still defined by Eq. (gj), but with given by Eq. (Q). It 
has to be pointed out that the interpretation of Eqs. (^) and ( ^ ) as the charge conservation law Eq. ( ^5| ) remains 
valid even for Ilgp 7^ 0. The only difference is that the expressions ( ^6|) for current and density 



are now more complicated and contain Ilgg . Nevertheless substitution of Eq. (^5|) in Eq. (^5|) indeed results in Eq. ( |34| ) 
with Mg given by Eq. Q). Hence, even in the most general case Hg3 7^ the position of zeros of Mg^ (poles of Mg) 
is gauge independent. 



A. Equation for the CG mode, its physical interpretation and gauge independence 

It can be checked that in the limit K ^ and fl/ K 00 one of the collective modes is the plasma mode considered 
above. We are, however, more interested whether a sound- like mode {n/ K = const for K ^ 0) which would be similar 
to the BA mode in neutral superconductors can exist in the charged system. It is seen from Eq. (^) that if the ratio 
fl/K is fixed and one is interested in the low energy excitations with ^ 0, only the last term of Eq. (^) is relevant 
because n33(i4r) and ng3(i^) are ~ V~^{K.) in this limit. In this case the equation A^^^(ri, K) = after the analytical 
continuation if2„ ^ ft + iO reduces to 



A'^^in, K)K^K, = I A«/'(f7, K) - ng3(».K)ng (|7,K) \ ^^^^ ^ ^^g^ 
y ll33(ii, K) J 

so that the detailed form of the Coulomb interaction becomes irrelevant, as was mentioned above. The solutions of 
Eq. (Eq) are gauge independent because for K — > this equation is equivalent to the gauge independent equation 

MgH^.K) = 0. 

It is interesting that practically the same arguments about the gauge independence of the CG mode studied using 
the formalism of Refs. |o[|ii], can be made. The corresponding equation for the CG mode is 

n22(f^,K)-H22(0,O) = 0, (47) 

where 

n22(f^, K) = n22(f^, K) + n23(f^, ^) ^^yj^^^^^^^ K) "^^^"' ^^^^ 

H23 is given by Eq. ( p8[ ) and the definition of 022 is given in (in addition to two obvious T2 matrices it also 
contains the factor 7^(k) because d-wave pairing is considered). One can show that in the limit K ^ equation ( |4^ ) 
reduces to the Ward identity 

h2/l^ + n22(f^,K^0)]H33(f^,K^0)+n23(f7,K^0) = 0. (49) 

Exactly this identity was recently proven in [2^ (to be precise, s-wave pairing was considered in |p^ ) using the gauge 
independence (charge conservation) arguments (see e.g. p7| ) and the mean field gap equation 

| + n22(o,o) = o. (50) 

As was already mentioned the last equation has to be explicitly used in the formalism of |l^, |ll|, |l^ . 
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It is possible to establish a link between Eq. (|4^) and a simple and transparent interpretation of the CG mode 
suggested by Schmid and Schon [38| (see also Chap. 13 in [^). Comparing the expression for current in Eq. ( ^ ) and 
Eq. (H) one can notice that since IIq^{K) disappears in the limit K — > 0, Eq. ( |46| ) is just a condition 

V • j - 0, j„(*r!„, K) = ^A"''{K)tKpe{K) (51) 

which implies that if there is a supercurrent in the system, it should be compensated by some normal current 
minimizing the total current. Exactly the same condition which relates the CG mode to a counterflow of supercurrent 
and normal current is discussed in [||. Thus using the two fluid picture one may say that A"^{K) in Eq. ( |46| ) consists 
of the supercurrent A"^ and normal A"'' (A') — A"^{K) parts, respectively. 

Although this counterflow resembles second sound in He'', it was pointed out in the earliest studies (see also 
^) that the CG mode is not the second sound which is a hydrodynamic mode since in the CG mode the normal 
fluid and superfluid are not in the local thermodynamic equilibrium and this mode is not a hydrodynamic mode. 

It has to be stressed that even though the solutions of Eqs. ( ^ ) and (^) are gauge independent for K ^ 0, the 
existence of such solutions is not required by the gauge invariance because no general statement can be made about 
the polarizations II33, A°'^ and Ilpg for arbitrary values of Vl and K. Thus, in contrast to the BA and plasma modes 
the existence of which is guaranteed by the Goldstone theorem and Anderson-Higgs mechanism, the CG mode 
does not obey any theorem and its existence is a fortunate result of many subtle features of the system dynamics. 



B. Velocity of the CG mode 



It is difflcult to solve Eq. ( ^ ) analytically for d-wave pairing and even numerically a more simple equation for real 
part, ReA"^(r2,K) = is usually considered ||ll|, |l^. This signiflcantly simplifies its solution, but in general this 
approximation can be justified only a posteriory , when the imaginary part is estimated. It is possible to study this 
equation in two ways. The first way is to extract the dispersion law r2(K) for, in general, arbitrary K. The second 
way is to find the velocity of the CG mode, vcg{4>) = in the limit K — > 0. The extraction of the dispersion 

law is more sensitive to the approximations whic h w er e ma de in the calculation of the polarization operators. In 
particular, usage of the approximated expressions ([A^), (|A4[ ) and (A5), which neglec t th e pair brea king for > 2Ao, 
introduces the restriction vpK < Aq. Moreover, for the polarizations Eqs. (C4), (C7) and (Cll) calculated using 
the nodal approximation (see Appendix ^) the condition of smallness of K becomes even more strict, so that here 
we study only the equation for vcg- As also discussed in Appendix ^ the nodal approximation is valid for T <^ Aq. 
This is the reason why in the present paper only the temperatures T < 0.6Tc are considered, where Tc is defined by 
Eq. ( C1C| ). Nevertheless, the presented formalism allows in principle to study the phase fluctuation structure factor 
up to Tc if no additional approximations are made in calculation of the polarizations. 

Thus for the illustrative purpose and comparison with p2| we also solve numerically the equation 



ReA 



Q/3 



n 

K 



Kr.K. 







(52) 



(it is rewritten in dimensionless form (Dl) for the numerical work) instead of Eq. (^) and the results are presented 
in Fig. ||. We stress that for — 7r/16 we obtain two solutions: vco ~ 0.83wf and vcg ~ 0.56wf- These results 
are in fact in a very good agreement with the results shown in Figs. 2 (a) and (b) of |T^. In particular, we also 
obtain that the velocity of the CG mode is practically temperature independent and its value is well described by 
VCG = vp cos((/> — 7r/4) (or v'(jq = vp cos((/i — 7r/4 + 7r/2) for the other pair of nodes which is rotated with respect to 
the first pair by 7r/2). Indeed, this equation gives 



VCG = VF y. 



0.707 (0.711 


at 


T = 


0.5Tc) : 


</) = 


(antinodal direction) 


0.831 (0.833 


at 


T = 


0.5T,) : 


(j) = 7r/16 




0.555 (0.562 


at 


T = 


0.5r,) : 


(j) = 7r/16 


(other pair of nodes) 


1 (1.00 


at 


T = 


0.5r,) : 


(j) = 7r/4 


(nodal direction) 



(53) 



Interestingly we obtain that the CG mode disappears even at somewhat lower temperature than in |12[] , T = O.lTc. 
This lowest value of T when the CG mode still exists is, however, more the result of numerical solution of Eq. (^2|) than 
a real threshold temperature, because the peaks of the structure factor considered below disappear rather gradually. 



C. Structure factor 



The knowledge of the structure factor S{Vl, K) is even more important than the value of the velocity, vcG^ because 
this factor contains the information about the damping of the CG mode [Bol. Furthermore, as we already mentioned. 
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this is the quantity which is measured in the Carlson-Goldman experiment . 

One of the main results of the present paper is that we obtain a closed analytical expression for the structure factor 
(3^ ) for the clean d-wave superconductor which is just a substitution of 1133, 1X03 and A given by Eqs. (|C^), ( Cll ), 
( CTI) in Eqs. (Q) and (Q) (its inverse imaginary part) in Eq. (|3^), respectively. The dimensionless form of Mg^ which 
convenient for numerical evaluation is given in Appendix ^ (Eqs. (Q and (Q). In our formalism Me(ri,K) for 
a fix ed va lue of vpK depends on the following dimensionless ratios: (i) ao — vf/va, (m) T/ep = T/Tc{Tc/ef) (see 
Eq. ( C10| ) where Tc is expressed via ep), {Hi) VpK'^/ujp = VpR'^/ej^^ej^/cUp) and the angle ((> which characterizes the 
direction of K with respect to the Fermi surface (see Appendix^. The ratio VpK'^ /ujp which was absent in Eq. ( |5^ ) 
is now present because the full expression for M^^ contains the Coulomb potential Vc(K) (see its representation 
in terms of vpK and tOp in Eq. D4) which we chose in the simplest 3D form. It is indeed easy to see that since 
vpK <C LOp, the detailed form of is not important and here we compute S{fl, K) including the Coulomb potential 
only to demonstrate this explicitly. 

Despite the gauge dependent definition of the pair susceptibility Eq. (||) , the peaks of the structure factor Eq. ( |3^ ) 
considered below which are associated with the singularities of M0(fl,K.) (zeros of A4g^{Q,'K)) can be regarded as 
gauge independent in the sense that the position of these singularities is gauge independent as was argued above. 
Furthermore, as shown in |^5| a more general expression for Aie derived from Eq. (2^) does depend on the gauge 
fixing parameter A, but in a such way that the position of the pole of Mg and its residue are gauge invariant. This 
justifies the use of the structure factor S'(r2,K) which is expressed via ImA^e(ri,K). 

In Fig. 1^ (a) - (c) we show the structure factor calculated using the analytical expressions mentioned above for 
different temperatures. The position of the peaks defined by the ratio ri/ypK is well fitted by Eq. (p3|). In particular, 
for (f) — 7r/16 two peaks are seen. This allows to associate the origin of these peaks with the gap nodes. As expected 
the peaks are getting less sharp and higher as the temperature increases. The width of the peaks also depends on 
the direction of K: for a larger value of </) (0 < < 7r/4) the corresponding peak is wider. All these results are in 
agreement with Fig. 3 from |12| , but due to the analytical character of the calculation the subtle peak features have 
a better resolution. 

To argue that these peaks in Fig. ^ are indeed due to the density-current coupling in Fig. ^ we show for comparison 
the structure factor which was calculated setting IlgQ = 0. The disappearance of the peaks confirms our claim that 
the CG mode demands nonzero density-current coupling. 

This procedure of setting Ilgg equal to zero is in fact very convenient in clarifying the origin of the peaks, because 
even for ng3 some peaks can be seen because both 1133 and IIqi^ have the square root singularities discussed in Sec. |rv| . 
For example, for 7r/16 < < 7r/4 (see Fig. ||) the lower peak becomes even sharper, but in contrast to to Fig. ^t 
does not disappear when we set Ilpg = 0. 



D. Manifestations of the CG mode in the classical response functions 

The structure factor considered above is so important for the investigation of the CG mode because for supercon- 
ductors there is no a classical laboratory field that couples, to and in the static limit is the thermodynamic conjugate 
of the order parameter. The reason that there is no laboratory conjugate field in the superconductor and superfiuid 
cases is that these order parameters are off-diagonal in number space Q . Nevertheless it is interesting to investigate 
whether the CG mode can manifest itself in the "classical" correlators when one goes beyond the static limit. 

Let us firstly consider the form Mt^ for the electric potential. Substituting the elements of the matrix Eq. (|2^) in 
Eq. ( ^ ) we arrive at the following expression (compare with Eq. (|37|)) 



One can verify that both plasma and CG modes are present in the equation A4^^ — 0. 



(54) 



Although the equation Ai^^^ = is gauge independent, Ai^^^ itself as well as A4f, ^ is gauge dependent. This time. 



however, we can take instead of Ai^ truly gauge independent correlators Eq. (pQI). Note that the definition Eq. (|3£ 
with the source fields (Aq, A) clearly shows that the "classical laboratory field" is used to probe the corresponding 
response. Then repeating the calculation of Eqs. (^ij) and ( ^2|) with nonzero Ilgg we arrive at the gauge independent 
density-density 
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and current-current 

correlators. 

In Fig. H we show the spectral density 

B(f],K) = llnix(r!,K) (57) 

TT 

calculated for the density-density correlator Eq. (^). The plasma mode (see Fig. ^ (b)) and the CG mode (see Fig. ^ 
(a)) as proves Fig. ^ (c) are clearly seen in the density-density correlator. We note that for (/) = 7r/16 the lower peak 
does not change strongly showing that it is due to the density-density correlations and not the CG mode. 

As the value of K decreases the relative weight of the CG mode with respect to the plasma mode becomes smaller 
and smaller as shown in Fig. |^. Finally in the limit 

xi^^n,K) Ao^PK^Kp ' ]^ = const, K ^ 0, (58) 

the CG mode disappears from the density-density correlator. This agrees with the statement of that the density- 
density correlator in the limit K — > is dominated by the plasma oscillation only. 
Treating the current-current correlator Eq. ( |56| ) in the same way one obtains 

X'^'mn, K) ^ (a"/' - ^^J§^) , § - const, n ^ 0. (59) 

Taking into account that the structure of A"'' is A"^ — ASap — FvFaVpp one can check that the term with F cancels 
out from the transverse correlator. This shows that the longitudinal CG mode cannot be seen in the transverse 
current-current correlator and its presence does not have any influence on the Meissner effect in the limit K ^ 0. 

VII. DISCUSSION 

In the present paper we have shown that the 2D model of clean d-wave superconductor predicts the existence of 
the CG mode in a wide temperature region down to T = O.lTc. This is done using the analytical expression for the 
structure factor which has peaks associated with the CG mode and solving numerically the equation for the CG mode 
velocity, vcg- All our results are in a good agreement with the paper h% where a similar model has been studied 



numerically using the formalism of It was also shown in |12 that in contrast to s-wave superconductors where the 
impurities supressing the Landau damping result in more favourable conditions for the observation of the CG mode, 
in d-wave superconductors the CG mode disappears in the dirty system. 

Thus the main physical question is whether our prediction of the CG mode in a clean d-wave superconductor 
is relevant for HTSC cuprates which are very complex compounds. Recent measurements [^T| done in high-purity 
YBa2Cu307 crystals show that the in-plane mean free path I increases to ~ 1/xm below T = 2{)K{k, Q.22Tc). This 
suggests that these systems are deeply in the clean limit I ^ ^ lOA (^ is the in-plane superconducting coherence 
length) and the model we considered can be applied to describe the CG mode in these compounds. Using the value 
of the Fermi velocity vp ~ 1-8 x lO^cm/s from we predict that the velocity of the CG mode is expected to be 
within the range vcg — (l/%/2 — ^)vf ~ (1-3 — 1.8) x lO^m/c depending on the direction of K which is one or two 
orders of magnitude faster than the velocity of the CG mode observed in conventional superconductors. 

From the theoretical point of view there are many questions which deserve further theoretical investigation. First 
of all, it would still be important to consider the influence of impurities and inelastic scattering by antiferromagnetic 
spin fluctuations within the proposed formalism. It is also interesting to estimate the widening of the structure factor 
peaks which is discussed above Eq. (|l8|). This widening may become important if temperatures T < Tc are considered. 

The discovery of the CG mode in conventional superconductors led to much deeper understanding of superconduc- 
tivity, so we hope that the investigation of the same problem in HTSC would also increase understanding of these 
complex systems. 
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APPENDIX A: GENERAL EXPRESSIONS FOR POLARIZATIONS Hij 

The general expressions for the polarization functions 11^ are (see e.g. [^l], p^ ) 



(27r)2 \ 2 V E-E+ 

1 / g_g+±A_A4 

2 \ E^E. 



1 



E+ + E- 



[l-nF{E^)'nF{E+)] 



E, 



E^ + if2.„ 



E, 



E^ - ia 



E++ E- - iQn _ 

[nF{E_)-nF{E+)]\V±(k), 



(Al) 



and 



y±(k) 
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J (27r)2 








\2E. 


2E^ 



i'FQ(k)i;F/3(k), 
1, 



2E, 



2E^ 



1 



1 



E++E^+ iQn E++E^- in 
1 



[l~nF{E^)-nF{E+)] 



(A2) 



E+- E^+ ifln 



nFiE^)-nF{E+)]\vFc{^) 



where C± = ^(k±K/2), E± = £:(k±K/2) and A± = A(k±K/2). One can also check that ri^(il7„ , K) = n^3(if]„,K). 

The first and second terms in Eqs. (Al) and (A2) have a clear physical interpretation p^. The first term propor- 
tional l — nF{E^) — nF{E+) gives the contribution from "superfluid" electrons. The second term gives the contribution 
of the thermally excited quasiparticles (i.e. "normal" fluid com pon ent). This is the term responsible for the appear- 
ance of the Landau terms in the effective action (see e.g. 18, |2j]). The imaginary part of these terms is the only 
source of damping of the phase excitations in the clean system considered here. The physical origin of this damping 
is due to the scattering of the thermally excited quasiparticles from the phase excitations. 

Since in w h at fo llows we are interested in the limit 11, vfK <C Aq we may safely neglect fl and K in the first terms 
of Eqs. (^, (|A2|) and write 



n33(f^,K) 



(27r)2 [E^ik) 



AHk)^^^^E{k) , ei^) 2{E+-E_ 



:[nFiE^)-nFiE+)] 



n;?o^(i7,K) 



d^k 



(27r)2 {E+-E_)^-n^ 
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2T £;2(k) {E+- E^Y -n?' 

UFiE-) - nF{E+)]vFaO<-)vFi3 0<.), 



2{E+ - E. 



(27r)2 E{k) {E+~ E-f -9? 
This approximation is in agreement with the one used in Appendix C of 111 



[nF{E^) - nF{E+)]vFc.{\^). 



(A3) 



(A4) 



(A5) 
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APPENDIX B: NODAL APPROXIMATION 



The key values which are necessary for evaluation of the polarizations Eqs. (A3), (A4) and (A5) are the differences 
— E_ and np{E^) — np{E_). Expanding in K, 



E(k + K/2) - E(k - K/2) = Vg(k)K , 
where the quasiparticle group velocity is given by 

1 



E{k) 



[e(k)vf + A(k)vA] 



(Bl) 



(B2) 



It is obvious that due to the gap k-dependence Eq. (Bl) differs from the s-wave case ||2J] by the second term To 
perform the calculation analytically it is useful to rewrite Eqs. (Bl) and (B2) in terms of the nodal approximation 
described in detail in (see also Sec. VII of |18| where the imaginary parts of these polarizations were calculated.) 
In particular, using this approximation one has e(k) ~ = vpki = pcosS, A(k) ~ p2 = VAk2 = psinS with 
E(k) ~ p = yjjf+jjf = \/vpkf + f^fc2i where the quasiparticle momentum k — (fci,fc2) is written in the nodal 
coordinate system ki, k2 associated with j-th node (j = 1, . . . ,4). (Note that the angle S was denoted in |jl^, Q as 
if). Then 



VgK = vpKi cos(5 + VAK2smS = Fcos(5 - -i/;) , £;(k± K/2) < Aq . 



(B3) 



where the momentum K — {Ki,K2) of 6'-particle is also expressed in the nodal coordinate system ki, k2, so that 
Pi EE vpKi = PcostP, P2 ^ y^K2 = PsinV' and P = ^/{P^Y + {P^f = y/vj,K'( + viK^ (We denoted the 
components of P as P^,P^ to make them different from the node label Pj used in what follows.) Finally, we can 
approximate the difference npiE^) — npiE^) as 



nF{E+) - npiE^) 



dnpiE) 
dE 



dnpiE) 
dE 



Pcos((5 — ip). 



(B4) 



Recall also that in the nodal approximation the integral over the original Brillouin zone is replaced by the integration 
over 4 nodal sub-zones as 



(27r)2 



4 

E 



dkidk2 
72^ 



4 

E 



d^p 



(27r)^Ui?WA 



E 



, = 1^0 



pdp 
2TrvFVA 



dip 
2^ 



y/TTUpUA/a. 



(B5) 



It is necessary to underline that the nodal approximation is designed for the low temperature regime T <C Aq. Thus 
all polarizations derived in Appendix ^ are applicable for T <C Aq. 



APPENDIX C: CALCULATION OF Hij AND EQUATION FOR n 



After the substitution Eqs. (p33|), ( |B4[ ) in Eq. (A3) and integration overp (see |1^) Has is expressed via the integral 

'■2'^ dS cos^Scos'^{S~iP) 



I 



27r cos2(,5- V)-62 



(CI) 



with b = |il|/Pj. (Note that we omitted the sub-zone index j in tp and b.) This integral can be calculated using the 
table integral (3.682) from for b^ > 1 giving 



1 |!,|3 

/= -+b^cos2il; ' ' cos^7A + |b|v/b2-lsin^V>. 

2 VlF^ 



(C2) 



Then using the analytical continuation in the region b'^ < 1, we finally arrive at the result 



ln2 T 1 



TT VpVA 1 2 P, 



- -f -2 cos2?Aj - 



—pj cos tpj — —pj ^ 1 sii^ V'. 



+ i- 



^2 



cos + 1 - 2 

3 \ 3 y 



sin^ ipj 



8 1 



\n\ 



(C3) 
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where k{T) denotes the term which originates from the first term in the braces of Eq. ( |A3| ) which cannot be accurately 
calculated using the nodal approximation [l^ . It is easy to show, however, that K = m/TTa.tT — for the continuum 
s-wave pairing model. This can be used to write 1133 in the dimensionless form: 



TTl ~ 

n33(fi,K) = -n33(o,K), 



where 



(C4) 



n33(f^,K) 



^ln2T 



(C5) 



where • • • are the terms in the braces in Eq. (C3), ep — mvp/2 is the Fermi energy, ajj = vf/va (see Sec. II) and 
we assume that k 1. One can also check that the imaginary part of 1133 coincides with the expression calculated in 
1^ integrating S{Pcos{S - ip) - il). 
In the same way we arrive at the expression for the current-current polarization function Eq. (|A4|): 



^ In 2 T Pj cos^ 



TT VpVA 



1 



02 



- 1 



-If 



1 



p. 



(C6) 



The zero order superfluid stiffness Aq^ = Sapn/m for the continuum (mQ/3(k) = mSafj) system at T = 0, where n 
is the total carrier density which, of course, coincides with the density of the neutralizing background. Using the 
expression ep = Tm/m which is, strictly speaking, valid only for the 2D systems with the quadratic dispersion law, 
we can also rewrite the superfluid stiffness A^^ in the dimensionless form 



A 



TT 



Ao 



4 

E 



T 

In 2 — an, 



p2 

^ COS^ t(jj 



(C7) 



As was already mentioned after Eq. (C3), the terms which contain the averaging over the Fermi surface (see e.g. 
Eqs. ( p^ ) and (|5|)) cannot be accurately calculated using the nodal approximation Q. Thus, in general Aq as well 
as ii should be considered as a free parameter of the model. In particular, decreasing the value of Aq it is possible to 
describe a lowering of the zero temperature superfluid stiffness in HTSC. Nevertheless, for the numerical computations 
we will assume that Aq '--^ 1. It is easy to obtain (see e.g. that the static, zero momentum bare superfluid 

stiffness 



A = — Ao - 21n2an — 

TT V €F 



(C8) 



and that the velocity of the BA mode 



A 



vp Aq — 2ln2ap)T/ep 



033(0, K) V2' 



(C9) 



so that for Aq = k = 1 the BA mode velocity v{T ^ 0) = vp/\/2. Using Eq. (C8) one can estimate the temperature 
of Berezinskii-Kosterlitz-Thouless transition from the equation Tc — 7r/2A(Tc), which gives 



epAo 



2(l + ln2a£,)' 



We will use this definition of to express the te mpe rature T in the units of and ep. 
Finally, we obtain that the expression for Eq. (A5) 
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Pn 



(Cll) 



where we put inside the braces the dimensionless part, HQ^Ka/ K 
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APPENDIX D: EQUATIONS FOR vcg, STRUCTURE FACTOR AND TRANSFORMATION TO THE 

GLOBAL COORDINATE SYSTEM 



Substituting Eqs. (C4), ( |C7| ) and ( |Cll| ) in Eq. (^6|), we obtain the equation for CG mode in the dimensionless form 
which is convenient for numerical investigation 



A = 



2e% n. 



33 



The whole expression (tt4h for Alg can also be written as 



where 



n33(f^,K)- 



n 



-ff" (n K]— + llaP^^^iEl 



(Dl) 



(D2) 



(D3) 



where the dimensionless polarization functions II33, and A"'' were made from the full polarization functions 11, 



33 J 



11^3, A"!^ (see Eqs. (pi) and (|1)) in the same way as the polarizations Eqs. (|C4|), (|C7|) and (\Cll\). The only difference 



is that these full polarizations include the Coulomb potential V"c(K) (for simplicity we take the 3D potential) which 
for our purposes is convenient to rewrite as 



2n ujI 



(D4) 



where ujp is the plasma frequency defined after Eq. (pq). 

Although the local nodal coordinate systems [Pj , ipj) are very convenient for calculating of the polarization functions 
Eqs. (p5D, (C7) and (Cll), the final expressions for them and, for example, Eq. (D2) have to be calculated in the 
global or laboratory coordinate system (K,(f>). It is convenient to measure the angle </> from the vector k^, so that 
(f> — corresponds to the corner of the Fermi surface (see e.g. Fig. 1 in ||l^) and the first node is at = 7r/4. Thus 
the transformations from the global coordinate system into the local system related to the j-th node are 



P 



. K^jvl cos2 (0 - J + |(^- - 1)) + vl sin2 (0 - J + \{3 - 1)) , 

VpK f IT IT \ V^K ( TT TT \ 

COS = COS - 4 + 2 ^ ' '^^"'/■'J " ~p~ V"^ ~ 4 + 2 ^ "^V ' 



(D5) 
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FIG. 1: Real part of n33(fi,K) (or more precisely n33(f2,K) = (7r/m)n33(n, K), see Eq. (|C^)). We put T = 0.5Tc (T^ and the 
angle are defined in Appendices ^ andj^, respectively) and an ~ 10. 
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FIG. 2: Temperature dependence of the velocities of the BA mode (dashed line - 1), the CG mode for (j> — (full line - 2) and 
4> = 7r/16 (dot-dashed lines - 3). We put an = 10 and the velocities are expressed in units of vf- 
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FIG. 3: Structure factor S{n,K) in arbitrary units: (a) T = 0.2Tc, (b) T = OATc and (c) T = O.GTc. We put ujp = 0.5eF, 
vfK = O.OleF (i.e. vpK = 0.02^;^) and an = 10, so that Tc ~ 0.063eF. 
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FIG. 4: Structure factor S{fl, K) in arbitrary units for T = OATc- The rest of the parameters is the same as in Fig. ^ except 
for the fact that the density-current coupling is "switched off" {TIqs ~ 0.) 
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FIG. 5: Structure factor S{Q, K) in arbitrary units for T — OATc- (a) Iloa 7^ 0, (b) Iloa — 0. The rest of the parameters is the 
same as in Fig. pi 
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FIG. 6; Spectral density B(r2,K) in arbitrary units for T = O.STc. (a) < Q./vfK < 1 and / 0, (b) 49 < Q./vfK < 51 
and / 0, (c) < Q/vfK < 1 and Tl^s = 0. We put ujp = O.Sef, ufA' = O.OleF (i.e. ujp = 50i'f-K') and an = 10. To display 
the (5-like plasma peak seen in Fig. ^ (b) a small imaginary part is added to uip. 
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FIG. 7: Spectral density B(Q., K) in arbitrary but the same as in Fig. gunits for T = O.STc We put 
(i.e. = SOOfF-ft^) and qd = 10. 



